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Abstract 

We discuss Brownian motion from the more elementary viewpoint 
presented by Einstein in 1908. Later developments and applications 
are briefly reviewed. 

One of the celebrated papers of 1905 deals with a quantitative study of the 
random motion. A couple of years later Einstein wrote a simplified account of 
the seminal work and we will present this version in this article. A direct test 
of the molecular kinetic theory was his obsession during his student days. The 
analysis of motion of suspended particles in a liquid was an attempt to make 
observable predictions that when tested would establish, beyond doubt, the 
correctness of the molecular point of view. 

The analysis begins by considering a non-uniform solution in a container 

with a semi-permeable membrane of thickness Ax in the centre, as shown in 

Fig 1. We need to consider the force balance in the region of width Ax. 

The pressure will be different on either side of the membrane because of 

dP 

the differing concentration. This will cause a force ———A/S.X to the left, 

ox 

where A is the cross-sectional area. Clearly, if the pressure decreases as x 
increases, then the direction of the force will be to the right. If the local 
density gradient is negative, that is the concentration is higher to the left, 
then there will be a drift of the solute molecules to the right. If the drift speed 
V and the shear viscosity of the solvent is rj, then there will be a resistive 
force Gn-qrv (Stokes Law) on the solute molecule, if it is taken to be a sphere 
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Figure 1: 



of radius r. These are the only forces in the horizontal direction and a steady 
state would require that the net external force is zero. 

If n{x) is the local density of solute molecules, then the total number in the 
volume AAx is n{x)AAx and the total viscous force will be Giri] r v n{x)AAx. 
The force balance then yields 

dP 

6nrirvn(x)AAx = ———AAx 

ox 



or 

1 dP 



67r?7 rn dx 
The solution pressure can be written as 



(1) 



P = nkT (2) 

at temperature T and thus 

1 dn 



Girrirn dx 



(3) 



The current of the solute molecules is the number of them crossing an unit 
area per unit time and hence can be written as 

j = nv (4) 

From Eq.Q, we find 

ksT dn 

^ = ~a IT 

DVT ?7 r ax 

The diffusion coefficient D is defined by the relation 
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Figure 2: 



and thus from Eq.® and Eq.® we have 

67r?7r ^^"^ 

We now look at the statistical property of the current. In the container 
shown in Fig.(l), we consider any cross section (the solid line in Fig. 2) 
AB and would like to calculate the current across the section due to the 
random hopping of the solute molecules. This random motion is produced 
by the random collisions of the solute molecules with the molecules of the 
solvent. In a time r, we assume the root mean square displacement of a 
solute particle to A. Consequently, a current across the section AB will be 
set up by those solute particles on the left of it that are moving to the right 
and solute particles to the right of AB that are moving to the left. Since 
the probabilities of a left or right jump are equal, half the particles within a 
distance A to the left of AB will cross AB in time r and half the particles 
within a distance A to the right of AB will cross it in time r. If rii is the 
density of the particles to the left and rir the density of the particles to the 
right of AB, then the number n of particle crossing an unit area of AB in 
time r is given by 

n = ^{niA-rirA) 

Hence, the number crossing unit area per unit time (which is the current) is 

1 A2 dn 

Comparing with Eq.(6), we have the second important conclusion of Einstein 

A^ = 2Dt (10) 
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Between Eq.d?)) and Eq.(ITU)), one has enough predictive power to find the 
diameter of suspended particles or Boltzmann's constant ks- 

The verification of Einstein's predictions was primarily the work of Perrin 
and his students. Early indications of the correctness of Eq.((Zj) and Eq. (|Tn|) 
came from the observations of Seddig, who took two photographs of an aque- 
ous suspension of cinnabar on the same plate at an interval of 0.1 second and 
measured the distance of corresponding images on the plate. He found that 
the distances at different temperatures were inversely proportional to viscos- 
ity as predicted. Perrin and his students followed the movements of single 
particles of gamboge or mastic under a microscope and recorded their po- 
sitions at equidistant time intervals by means of an indicating apparatus. 

Since the particle size was known, these observations yielded kB and since R 

R 

was known, one could get N = - — . 

ks 

Perrin explicitly established that the suspended particles were in thermal 
equilibrium with the solvent by studying the distribution of particles in a 
vertical column under the action of gravity. The total energy of a suspended 

particle at a height z above the base is h mgz and hence the number of 

particles between z and z + dz is given by 

/d^p p^ Tfigz 

where C is a constant and A is the area of the base of the container. The 
total number of particles is found from 

N^CAjileM-^)jeM-^) (12) 

From Eq. ljllj) and Eq. ()12|l we have 

This distribution was confirmed in Perrin's experiment (Fig. 3a). The 
displacement of individual particles had the typical from shown in Fig 3b. 
The path of the particle is an example of a fractal, a curve for which any 
small section resembles the curve as a whole. 

A different way of looking at the problem was derived by Paul Langevin, 
who was a friend of Einstein. Consider a molecule of mass M colliding 
elastically with another molecule of mass m. If the velocities of M before 
and after the collision are V and V' respectively and those of m are v and v' 
then 

MV + mv = MV' + mv' (14) 

and 

-MV^ + -mv^ = -MV'^ + -mv'^ (15) 
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Figure 3: 



From the above equations we immediately find that 

V = + (16) 

M+m M+m ^ ' 

We now imagine that the molecule M is being hit randomly by the molecules 
m from all directions. In each collision, the change in momentum is 

^ ' M + m M + m ^ ' 

The force on M due to the random collisions, is then clearly composed of 
two parts : 

i) a resistive part which is proportional to V . 

ii) a random part which is proportional to the velocity of the irregular mo- 
tion of the molecule m. The average value of the force is zero and the mean 
square value is proportional to fc^T, where T is the temperature because the 
mean square velocity of the gas molecule is proportional to ksT . 

The equation of motion for the molecule can be written as 

dV - - 

- = -n. + / (18) 

where F is the relaxation rate and / is a random force with 

</, > = (19) 
</*(t)/,(0> = 2a<5(t - (20) 
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The above arguments tell us that a oc ksT and if M is a spherical solute 
molecule moving through the fluid, then F = ^ where r is the radius. 

The relation between a and diffusion constant D is important for the estab- 
lishment of thermal equilibrium. 

The solution of Langevin equation can be written down as 

V,{t) = e^* f e^''m')dt' + V,(0)e-r* (21) 
Jo 

The mean square velocity is can be written as 

< V^{t) >=< e-^r* fe^^'mdt' fe^'"m")dt" > + 
Jo Jo 

V;2(0)e-2r* + 2 < r(0)e"^* / e^'' fi{t')dt' > (22) 

Jo 

The averaging is over the random force f{t) and using Eq. (jl9|) and Eq. (j2(Jj) 
we have 



< V^it) > = 2ae-2r* J e"^'' dt' + V\Q)e-^^' 



= f + mO)-f)e-''' (23) 

If the solute molecules come into thermal equilibrium with the solvent molecules, 
then they too will have a Maxwell Boltzmann velocity distribution in which 
case < vf >= and we see that ^ = The displacement Axj in time 
t can be written as 

Ax, = [ Vi{t')dt' 
Jo 

ft pt' 



f dt'e-^'' f citie^V4ti) + ^(l-e-^*) (24) 
JO Jo ^ 



This leads to 



<fAx,)2>=< C dt'e-^'' 







' f e^''Mh)dh [ dt"e'^'" [ e^''U{h)dh > (25) 
Jo Jo Jo 

[ dt'e-^'' [ e^''Mti)dti > 
Jo Jo 
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Using the moments of f{t) 



+ 



2a / dt' / rft"e-^(*'+*") 
Jo Jo 

/ dtiexp{2Tti)e{t' -t" 
Jo 



[ dtiexp{2Vti)e{t" -t') 
Jo 



r2 



■[I - exp{-Tt)Y 



dt' / dt" exp{-T{t" -t'))e{t' -t") 
Jo 

+ I dt' I dt"exp(-T(t' - t"))9(t' - t")] + ^-^[1 - expi-Tt)? 
Jo Jo J ^ 

[ dt' [ dt"exp{-T{t" -t'))+ [ dt' [ dt"exp{-T{t' -t 
Jo Jo Jo Jo 

[1 - exp{-Tt)f 



+ 



a 
f 

Y2 



')) 



a 



dt"[l - exp{-Tt")] + 



dt'[l- exp{-Vt')] 



+ 



Tt\2 



Y2 



(26) 



For extremely long times, the first term dominates and 

/A n9 2(7 

<(Ax,)2> = —t 

keTt 

M r 



GTirir 
2Dt 



t 



(27) 



where D 



Gurir 



is the Einstein relation as found in Eq.((7j) and Eq. ()10|) . 



There is yet another way of approaching the problem. This is through the 
Fokker-Planck equation which gives an equation for the dynamics of the 
probability distribution associated with the random process described by the 
Langevin description of Eq. ()18|) . If P{x,t) is the probability distribution as- 
sociated with the random variable x{t) which satisfies the Langevin equation 



dx ^dS ^, , 



(28) 



with the random force being Gaussian and having the correlation 

< f{ti)f{t2)>=2D6{t2-t,) (29) 
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Then the probabihty P{x,t) of the variable X{t) having a value x at time t 
satisfies the evolution equation 



In the case of a particle simply changing it's position because of random 

dS 

fluctuations we can drop the deterministic part of Eq.()28|) i.e — — and that 



implies a probability evolution equation 



dx 



If a; is xo at t = to, then the probability P{x,t, |xo,to) at a subsequent time 
is easily seen to be 

P(..,N.io)^ ^J^^_^^^ exp(-j|^) (32) 

We can derive the mean square displacement as 

<{l\xY> = < {x - xqY >= / dx ^ exp' ^ 



— oo 



^2nD{t-to) 2D{t-to)' 
= 2D{t - to) = 2DAt (33) 

which brings us back to Eq. ()10|) once more. 

This way of formulating the problem allows us to study an example of 
crossover behaviour in a straight forward fashion. Instead of allowing the 
randomly moving particle all of space to wander about , we restrict the 
motion in a region of length L (we will stick to one dimensional case, the 
generalisation is obvious). The particle reflects from the boundary walls 
placed at a: = and x = L. It is clear that the mean square displacement has 
to be of 0{L'^) independent of time as time goes on increasing indefinitely. 
On the other hand Einstein's relation tells us that for L — > oo the mean 
square displacement has to be order r. Thus, 

{/\xY ~ for r ^ oo at finite L. 



and 



{/\xf ^ r for L ^ oo (34) 

This is a classic crossover problem, the crossover being determined by the 
scales L and r. For » r it is the usual Einstein relation, while for 
<< r it is a size limited answer. The crossover function which interpolates 
between the two limits can be constructed by solving Eq.(jnH) in a finite 
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x=0 

Figure 4: 

geometry. The normalized version of Eq. ()31|) with the boundary condition 
dP 

that -— = (reflects from the wall, no current) at x = and x = L, can be 
ox 

written as 

Fix, t|0, 0) = - + - ^ cos(-^) exp( ) (35) 

n=l 

The mean square displacement is readily found to be 

< Axf > = + ^^^exp 36 

n=l 

L2 4L2 ^ 1 -n^TT^Dt, 4L2 ^ 1 y^^Dt^ 
= 5- / ^exp( — )H > ^exp( — - — ) 

n odd n eiieri 

(37) 

It is straight forward to check that 

< (AxY > — > 00 for t — >• 00 

and 



< {Axf > 2Dt for L^oo (38) 

To find the crossover function, we need to express Eq. ljHTj) in a closed form. 
This can be done by using the Euler-Maclaurin sum formula for the two sums 
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shown in Eq.(|37|). The sum formula is given by 



f{x)dx = h[^f{a) + f{a + h) + f{a + 2h) + --- + f{b - h) + ^f{b)] 



where Bo and Ba are the Bernoulh numbers with 



(39) 



1 1 

Bo = - and Ba = 

6 30 



(40) 



Keeping only exp{- 



L2 



from among the different exponential decays and 



working to 0(-^) the coefficient of exp{- 



be written as 



L2 



the crossover function can 



< (Aa;) >= — exp( — — )dn 

3 TT"^ Ji 77,^ 

' - o)(l + ^^)exp( 7^) 



L2 



L2 



(41) 



Going for a stronger approximation, we can write the crossover function as 



'2 r 



< {/\xY >-- 



i-a 



6)t. 



TT 



■Dt. 



L2 



exp( 



L2 



(42) 



The numerical simulation of the finite size diffusion and it's comparison with 
our crossover function are shown in Fig.5. 

The Langevin equation approach has been of tremendous importance in 
extending Einstein's picture of molecular motion to several systems over the 
last one hundred years. One of the prime areas in this respect has been the 
dynamics near second order phase transitions. For over three decades, the 
nineteen sixties, seventies and eighties, critical dynamics had been one of the 
frontier areas of research in condensed matter physics. A typical example 
would be the liquid gas critical point which was found by Andrews almost one 
hundred years ago. In his extensive experiment on carbon dioxide, Andrews 
found that at low temperatures as one increases the pressure on a dilute gas, 
the volume decreased and at a particular density the gas condensed to a liquid 
and over a range of density there were coexisting liquid and vapour phases. 

The pressure had to be increased enormously to cause a of density after 

the whole gas had liquified. However, about a critical temperature Tc, the 
gas could not be liquified by application of pressure. Consequently at Tc, the 
coexisting phase diagram disappears and one has a single phase region. The 
transition at Tq from a one phase system to a coexisting phases system is a 



10 




second order transition and is characterised by large scale fluctuations. As 
one lowers the temperature towards the critical point, droplets of the liquid 
phase are formed due to the fluctuations and the droplets become larger and 
larger as the critical point is approached. If the characteristic size of the 
droplet is ^, the C, becomes inflnitely big aX T = Tq since the droplets then 
become of macroscopic size and at Tc the liquid and vapour phases coexist. 
This divergence is in general characterised by an exponent u and one writes 

e = eo(T-7br (43) 

where is a system dependent constant and f is system independent con- 
stant. 

We now imagine a thermal gradient applied to this system. The droplets 
diffuse to carry the heat current and as the size of the droplets increase the 
process is expected to become more and more inefficient. This is what we 
will be able to demonstrate by using the arguments that we used in the early 
part of the article. A temperature gradient produces a pressure gradient and 
a force per droplet which can be written as 

The viscous forces on the droplet is dnrj^v, which yields 

V = (45) 
Gttt; { dx 
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The heat current can be written as 



J. = TCpV (46) 

where Cp is the constant pressure specific heat per unit volume. The heat 
current by definition is 

where A is the thermal conductivity. From Eq. ()46|) and Eq. ()47p we find 



Cp 6nr] ^ 



(48) 



where — is the coefficient of heat diffusion. As expected diffusion becomes 
Cp 

very small as T approaches Tq because of the divergence of the correlation 
length ^. On the other hand, the specific heat at constant pressure diverges 
as Tq is approached and the divergence goes Consequently, Eq. ljlHj) leads 
to the remarkable result that the thermal conductivity of a fiuid diverges 
near the critical point since 

A = ^ oce (49) 

This spectacular result was first obtained experimentally by Jan Sengers in 
the early sixties and the theoretical result established in the seventies. It 
is interesting to see how the simple analysis carried out by Einstein can be 
invoked to understand a path breaking result found more than fifty years 
later. 

We now turn to another aspect of diffusion which has seen significant 
alterations in the course of last one hundred years. If we drop a lump of 
sugar in a beaker of water and let molecular diffusion cause the sugar to 
dissolve in water, then it is a common experience that it takes a long time 
for the sugar to mix uniformly in the liquid. We can use Einstein's relation to 
estimate the time. If L is the typical dimension of the beaker, then according 
to Eq. pO|) the mixing time is given by 

10^ cm^ 
~ 2D ~ 2 X 10-2 cmVsec 
^ 1.4 hrs 

for a beaker with L ~ 10 cm and a diffusion coefficient of 10"^ cm^/sec (typ- 
ical of sugar in water). This is an awfully long time and that is why liquid 
is stirred in order to make the sugar mix uniformly in a short time. If v is 
the typical velocity, then the mixing is now 

-'m = ^ (51) 
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Clearly the ratio of the two mixing times is given by 



r„ _ L 2D _ 2u D 

- i 

where v is the momentum diffusivity that is z/ = — (kinematic viscosity) 
and S and Re are two dimensionless number known as Schmidt number and 
Reynolds number respectively. The Schmidt number S = — is the ratio of 

V 

the two diffusivities D and v and is a number of order unity. The Reynolds 

vL 

number Re = — requires an analysis of fluid flow. 

V 

The dynamics of an incompressible fluid is governed by Navier Stokes 
equation which reads 

^ + (t;. V)t; = -^ + i/V'w (53) 

The terms on the right hand side are the forces-a force coming from a pressure 
gradient and a force coming from the viscous drag. On the left hand side, 
we have total acceleration. The first term originates from an explicit time 
dependence of the velocity and the second from the fact that the flow brings 
in and takes out fluid from an elementary volume and thus causes change in 
momentum of the fluid in the volume considered. The second term in the 
left hand side is nonlinear and the source of all the difficulties in the solution 
of Navier-Stokes equation. This is the term which becomes important when 
the velocity is large and causes instability of laminar flows and eventually 
gives rise to turbulence. How important is the term with respect to the other 

terms? The estimate of this term (called the inertial force) is — where v is 

L 

typical velocity and L a typical length. The estimate of the viscous drag is 
— and thus 



Inertial force f ^ 

Re (54) 



Viscousdrag L vv 

Lv 



V 



Thus the Reynold's number of a flow is an indicator of how important 
the nonlinear term is in comparison to the linear term. Larger the Reynold's 
number, the greater the role of the nonlinear term. We now see from Eq. ()52|1 . 
that the mixing time is going to get smaller as the Reynold's number is in- 
creased and in the limit of Re <— oo (fully developed turbulence), the mixing 
time Tra in the presence of stirring seems to go to zero which is an indica- 
tion of the fact that in this limit, Eq. lfTnjl needs to be modified. Turbulent 
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diffusion differs very strongly from ordinary diffusion and is the mechanism 
behind the fast dissolution of sugar in a stirred beaker of water. The study 
of diffusive behaviour has drawn a lot of attention over the last couple of 
decades. 

The dissolution by stirring certainly requires us to study a forced Navier- 
Stokes equation. This can be written as 

dv 1 

— f + VfsdpVa = daP + l^V^Va + fa 

at p 

where /„ is a random function (specified by correlations only) that mod- 
els the stirring of the fluid. The fluid is incompressible which means that 
V ■ = 0. It is important to note the energy balance of the above equation. 
The total energy of the fluid can be written as 



^ J VaVa ' 



E = - I VaVa d r 



dE 
'dt 



= - y" VaVbetadpVa d^ - ^ J VaOaP d^V + U J VaV'^Va d^T + j faVa d^T 

= JPVadSa^O 

where we have assumed that the velocity vanishes in distant surfaces and 
have used daVa = in the first step. In a similar manner 

J VaVpdpVad^r ^ \j '^l3{vaVa)d^r 



d,3{v'Va) dh 



^ j v\f3 dSp 



= 

For the viscous dissipation, we have 



J VaV^Vad^r = -V j {df3Va){di3Va)d^r 



Thus, 

dE 
'dt 



-V J {dpVafd^r + J faVa dh 



'dE 

We can maintain a steady state in the system ( — = ) if we balance the 



dt 

viscous dissipation (first term on the right hand side) by the second term 



14 



which represents the rate at which energy is fed into the system. In this 
global analysis, the nonlinear term plays no role. Their contribution is in the 
carrying of the energy across the different length scales. The dissipation term 
which carries the highest derivatives in the problem is clearly most relevant at 
short length scales (dissipation effective at molecular scale), while the energy 
input (no derivatives) is operative at the highest length scales. Consequently 
the energy balance picture for fully developed turbulence that was arrived at 
Kolmogorov involves injecting energy into the system at large length scales 
at a constant rate e, transferring at from large scales to short scales at the 
constant rate e and dissipating it at the shortest (molecular) scales by viscous 
action at the same rate e. For this picture to be consistent, the highest and 
shortest length scales need to be well separated. Now, the scale, L,at which 
energy is injected can be estimated from 

e = 777 OC — 

T L 

where v is the typical velocity and T is a typical time scale over which the 
energy containing eddies of size L can turn in response to the drive. The 
dissipation scale is I and has to be constructed out of e and the kinematic 
viscosity v. Dimensional analysis shows that 



I 



1/3 



Consequently, 



L L fv^ 

1j 



I (^)i 



3 

i/4 



Thus for high Reynold's number, the scales are extremely well separated.we 
can see a particularly relevant feature of turbulent flow from this simple 
analysis. In the limit of Re tending to infinity, the scale / has to be extremely 
small and if we try to construct derivative of the velocity field, that derivative 

will become singular, which is exactly what is needed to keep v J d^r{daVpf' 

finite in the limit of z/ — * 0. 

We now note that equations of the form 



du 
dt 

2 



r 

have solutions which depend on the combination — . The solutions can 

change by a scale factor if r and t are scaled by appropriate factors. If r 
is scaled by A and t is scaled by A^, then the solutions of u can change 
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only by a power of A. Such solutions are known as scaling solutions. It 
should be apparent that Einstein's result for Brownian motion fEq. ()10|) is 
a manifestation of this scaling solution. If the diffusion law needs to be 
changed for the stirred fluid, then there has to be a different scaling solution. 
Accordingly, we explore the scaling solution of 

, A V72 

We have dropped the pressure term because taking a divergence of Eq.(|^. 
we note that ^ 

P 



and thus the term da — can be written as da 

P 



daVpdpVo 



and the dimension 



V2 

of this term is identical to that of VfjOpVa- Consequently, the pressure term 
cannot modify a dimensional analysis. 
We carry out the following scaling 



X = 


\x' 


t = 


\H' 


V = 





In the primed variables, we have 



or 



If we impose 



3 

then provided u is scale dependent and becomes 

We have the same Navier-Stokes equation in primed variables 

^""^ + v'd'v' = 



a 



The non-linear term has had two drastic effects 

(a) z is no longer 2 and hence Eq.()10|) will have to be modified. 

(b) the viscosity has become scale dependent-it is no longer the molecular 
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viscosity. We have generated a turbulence induced scale dependent viscosity. 



The scaling exponent z is however unknown yet. To explore that further 
we explore the effect of the scale transformation on the energy flow per unit 
time e. We find 

Since the rate of energy flow is constant, we require e = e' and that flxes 
2 

= 1, leading z = -. We now explore the effect of this on the concentration 
3 

dc -» 

evolution. The scale transformation takes tt + ('^' V)c = DV^c to, 

at 

F)r' 

+ A(t-i)(tr ■ V')c' = DA-^V'^c' 

or 

Fir' 

£^ + A(t+--i)(v' ■ V')c' = DVV 

where 

D' = \^-'D 

The advective terms in the diffusion equation has thus induced a scale de- 
pendent diffusion coefficient, where D has to scale as L^~^. In Eq.(ITU)), we 
now need the structure 

where Dq has the dimension of a diffusion constant and Iq is a suitable length 
scale. In terms of t, L ~ t^ and thus, 



{/\xf oc Do^z^t (X t^ = t^ 

^0 

where we have used z = ^ in the last step. Consequently, for turbulent 
diffusion 

(Ax)2 oc t^ 

which replaces Eq. fllOp of Einstein and ensures in a stirred fluid, the mixing 
takes place much faster than it would otherwise. 

We end this article by considering a particular development which has 
occurred over the last decade. It brought to light a feature of the diffusion 
process which had not been noted before. The diffusion equation for a con- 
centration n{x, t) is a conservation law for the total concentration. In the 
region of thickness Ax in flg.l, the amount of diffusing material is n(x, t)A/S.x, 
where A is the cross sectional area. In time At, the change in amount of ma- 
dfiix^ t\ 

terial is — — ^ — AAxAt. This change can only be brought about the solute 
ot 
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current (this is the conservation law) and is the difference in the current at 
the two surfaces bounding the region of width Ax. The net amount of ma- 
terial flowing into the region in time At because of the difference in current 
at the two surfaces is —^AxAAt. The conservation law then reads 

I = -f = (55) 

ot ax ox'^ 

on using Eq.®. We can without any loss of generality rescale 'x' and 't' 
to set -D = 1 and we can write the diffusion equation as 



On d'^n 



dt dx"^ 

This solution is well known and can be written as 



(56) 



n{x,t)= G{x-x',t)n{x',0)dx' (57) 



where ^ 

G{x — x' ,t) = exp —{x ~ x'Y / At (58) 

Clearly, the solution is characterised by a single length scale which grows 
as t^/^. What was discovered by Majumdar et al about ten years ago is that 
there is another non-trivial exponent associated with this apparently solved 
problem. 

To understand the origin of this new exponent, we consider a class of 
initial conditions where n{x, 0) is a Gaussian random variable with zero mean. 
The question asked is : what is the probability po{t) that the field n{x,t) at 
a particular point x has not changed its sign till time t. For large time, this 
probability has to decay to zero and this long time decay is charactrized by 
the law 

Poit) cc t-' (59) 

The exponent 6 is the new non-trivial exponent that we talked about and 
is called the persistence exponent. It also turns out that Pnitijh) - the 
probability that the field changes sign n times between ti and t2 (ti > ^2) - 
is given by 



Pn{tl,t2 



t2 



t. 



(60) 



The primary feature of the calculation is the point that the Gaussian process 
n{x, t) is a Gaussian stationary process in terms of the new variable T = Int. 
This is followed by the central assumption that the interval between two 
successive zeros of n{x,t) can be treated as independent. 

The distribution of n{x, 0) is taken to be Gaussian and white, i.e the two 
point correlation is 

<n{x,0)n{x',0)>=6{x-x') (61) 
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We can now write down the correlator as 



<n{x,0)n{x',0) > = < dx' G{x - x' ,h)(f){x' ,0) / dx" G{x - x" ,t2)(j)ix" ,0) > 



dx' dx" G{x - x\ ti)G{x - x", ta) < (j){x', 0)(j){x", 0) > 



dx' G{x — x', ti)G{x — x' , ta) 



where we have made the use of Eq. (jUT|) in arriving at the last step. Using 
the Green function given in Eq.()58|). 



< n(x, ti)n(x', ^2) > 



1 



1 



dx' exp 



(x — x') 



/\2 



20? (ti + ts)^ 



(63) 



We immediately see that 



< n{x, tiY > 
If we now define a variable t) 



1 



2^/27111 
n{x, t) 



correlation function 



\/n{x,t) 



22 (27r)4 t-i n{x,t), then the 



< X(ti)X(t2) >= 



2(27r)2(tit2)4 

2v/?(tl+t2)5 



4tit2 



(64) 



This is not yet stationary process. To construct a stationary process we need 
to use the new time variable T = In t. Now 



a(7'i,T2 



4exp(Ti +T2) 



T2\2 



[I 



=(T2-Tl)^ 



sech{-(T2-T0} 



(65) 



Thus one has arrived at a Gaussian stationary process, the anticipated form 



of Po{ti,h) is Poiti,t2) 



ti 



e ^^'^^ In the new time variable the 



calculation of 9 is the calculation of decay rate. In the process of calculation, 
one needs to make a very important assumption - the interval between zeros 
of X{T) are statistically independent. 



(62) 
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The correlation function a(T2 — Ti) will not give us the information re- 
garding the zeros of X{T). For that we need to know teh correlation of 
the sign of X{T). Accordingly we define a function (j{T) = sgnX{T) and 
consider the correlation 

A{T) =< a{0) a{T) > (66) 

The correlation function is determined entirely by the distribution P{T) of 
the intervals between zeros. The trick is to find the P{T) from A{T) and 
then Po{T) from P{T). The first step is knowing A{T) and this is possible 
for a Gaussian stationary process , because for such a process 

A{T) = -sin-^a{T) (67) 

TT 

Since a(T) has been calculated in Eq. (jU3j) we know A{T). Now if Pn{T) is 
the probability of there being n zeros between and T, then clearly 

oo 

MT) = J2(-irPn{T) (68) 

n=0 

Looking at the interval between and T, if there is a zero at Ti, then to find 
Pi{T), the probability of there being just one zero between and T then we 
nedd to introduce Q{T), the probability that an interval to the right or left 
of a zero contains no further zeros. Clearly 

PiiT) = £ dT, Q{T,)Q{T - Ti) (69) 

where < T > is the mean interval size. Now P{T) is the probability of finding 
an interval of length T between two successive zeros. So, if we are interested 
in the probability of there being two zeros between and Ti (the zeros being 
at Ti and T2 ), then 

P2{T) = r dT, r dT2 Q{T,)P{T2 - T,)Q{T - T^) (70) 

< ^ > Jo Jt 

if the intervals are all independent of each other. This is the "Independent 
Interval Approximation" and generalising to the case of n zeros, 

PniT) = — ^ r^/Ti r dT^... f dT^Q{n)P{T2 - T - 1)P{T, -T,).. 

...P{T^-T,,_,)Q{T-Tr, 
for n > 

The probability that an interval to the left or right of a zero does not contain 
a zero is unity if the interval is of unit length and hence 

Q{T) = 1 - / P{t)dt 
Jo 
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which leads to 

g'(T) = -P(T) (72) 

Using Laplace transforms, if P{s) and Q{s) are the transforms of P{T) and 
g(T), then 

P{s) = 1 - sQis) (73) 

and from Eq. (f7T|) 

Pn{s) = ^^[Qismp{s)r' (74) 

Eliminating Q{s) by using Eq. (f7Hj) we get 

pn{s) = - mnmr' for n > i (75) 

00 



For n = 0, we use ^^p„(t) = 1, which implies ^^p,i(s) = - to write 

n=0 n=0 ^ 

1 °° 11 

^^o(^) = - - T.p^is^ = - - ttt^^^ - ^^'^^ ^^^^ 

n=l 



The Laplace transform of Eq. ljUSj) gives A{s) = ^^(— l)"p„(s). We use 



Eq.(f7^ and Eq-lfTHj) to perform the sum and finally obtain 

F{s) = l + ^^.s[l-sA{s)] (78) 

The sequence now is as follows: From Eq. ()78|) . F{s) is obtained in terms 
of A{s) and from Eq.(f77j). P{s) is obtained in terms of A{s). We not that 
A{T) has been calculated in Eq. (jU7j) and Eq. ljU^ . Once P{s) is known po{s) 
and Pn{s) [n > 1] are known from Eq. ()76|l and Eq. (f75|) . The inverse Laplace 
transform will yield pn{T). 

It is however simpler to do the following. We note that Po{T) ~ e^^ and 
hence po{s) has a simple pole at s = —9. Looking at the structure of Eq.lfTmi 
and Eq. (f77j) . we note that a simple pole of po{s) is a simple zero of F{s). For 
the A{s) corresponding to Eq. (p7|) . 



1-— y dTe-'^sin-^[sech-]'A 



(79) 



This F{s) has only one zero and numerically it is found at s = —0.1203 
giving 

e = 0.1203 
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A numerical simulation of the diffusion equation yields 6 = 0.1207. Thus, 
we see how a new nontrivial and totally unexpected feature of the diffusion 
equation has been found almost a hundred years after the equation was first 
used. 
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